Thermodynamic stability of neutral Xe defects in diamond 
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Optically active defect centers in diamond are of considerable interest, and ab initio calculations 
have provided valuable insight into the physics of these systems. Candidate structures for the Xe 
center in diamond, for which little structural information is known, are modeled using Density 
Functional Theory. The relative thermodynamic stabilities were calculated for two likely structural 
arrangements. The split-vacancy structure is found to be the most stable for all temperatures up to 
1500K. A vibrational analysis was also carried out, predicting Raman- and IR-active modes which 
may aid in distinguishing between center structures. 
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I. INTRODUCTION 

There are over 500 optically active defects in 
diamond 1 - - — , few of which have had their structures 
completely determined. These defects' utilities range 
from potential bulk dopants to achieve superconductiv- 
ity in diamond 4 ^ to bulk emissive properties in light- 
emitting diodes 6 . Other example defects include the 
neutral- vacancy GR1 7 , nickel-based NIRIM-1 and -2&, 
boron acceptors 9 -, as well as the only room-temperature 
stable single-photon source defects such as the well- 
studied nitrogen- vacancy center (NV)i£, the nickel-based 
NE8 cente r 11 ' 12 , silicon- vacancy (SiV)^, the unidentified 
734 nm center— and the possibly chromium-based MHz 
emission center—. In addition, the nitrogen-vacancy 
(NV) center has also been shown to act as a possible 
qubit for quantum computing^ 6 -, as a magnetometer— ~—, 
and to be useful in probing quantum decoherenc o 21 ' 22 . 
It has also been used to demonstrate Quantum Key 
Distribution" using the BB84 scheme^. The optically- 
active Xe defect, the subject of this work, was only re- 
cently characterized 2 ^— and is of interest as a near-IR 
emitter with two zero phonon lines at 794 and 813 nm. 

Optically active defects in diamond have high poten- 
tial for technological applications, hence it is important 
to also consider their thermodynamic stability over the 
range of temperatures to which they are likely to be sub- 
jected, including during their creation. Information on 
the relative stability of defects is also paramount in opti- 
mizing their manufacturing conditions. Atomistic mod- 
eling of these types of defects offers several advantages: 
thermodynamic stabilities (both particular and relative) 
of defects can be evaluated, predictions of their vibra- 
tional spectrum made (which can aid in identifying the 
defect structure) and examination of potential electronic 
or optical properties may be undertaken. 

The observation in Refs i 27 ' 28 of a Xe-based defect in di- 
amond emitting at 794 and 813 nm makes the case for de- 
tailed ab initio studies of such centers even stronger. The 
site symmetry of that defect was established as trigonal 
by their polarization study. This defect is likely to be 
thermodynamically driven due to the experimental con- 



ditions under which it is formed (180-500 keV implanta- 
tion with 800-1400°C annealing), but little else has been 
published upon Xe-related defects in diamond other than 
recent studies indicating that Xe chemically reacts with 
the diamond lattice^, and that the defect is comprised 
of one Xe atom (with assorted vacancies) 28 . 

Density Functional Theory (DFT) 30 i 31 has provided 
important insight into the nature of such defect a 29 ' 32 " — 
and is excellent at predicting ground-state proper- 
ties. The stability of various morphologies 3 ^, mod- 
eling of dopant a 34 ' 36 ' 37 , and semiconductor band-gap 
modification 3 — have been characterised. 

As a material, diamond offers several other attrac- 
tive properties, including the largest optical transparency 
window, Young's modulus, and the advantage of a chemi- 
cally inert host lattice^ 9 .. The combination of this excep- 
tional host material with its plethora of dopants offers 
much scope for investigation and development. The lack 
of definitive structure models for many of these defect 
centers requires attention. 

In this paper, as the true geometric structure of the de- 
fect has not been identified by experiment, we use DFT to 
study the thermodynamic stability of neutral Xe- vacancy 
defects in diamond which both have the same symmetry 
as those observed by Bergman et a/. 27 ' 28 , and contain one 
Xe atom 2 ^. We specifically consider a substitutional Xe 
atom with either one (Xe s -V) or three (Xe s -3V) vacant 
lattice sites directly adjacent to the Xe atom. The paper 
is organised as follows: in Section |TT] we discuss defect 
formation free energies and how they may be compared; 
Sec. IIIII details the DFT parameters used in calculations; 
Sec. IIVI presents relaxed defect structures; in Sec. IVII 
vibrational results are discussed; the relative defect for- 
mation free energy is shown in Sec. [Vj and conclusions 
are drawn in Sec. IVIII 



II. DEFECT FORMATION ENERGY & 
THERMODYNAMIC STABILITY 

To compare formation free energies between chemical 
structures requires accounting for the changes one must 
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make to transform from one structure to the other. The 
defect formation free energy (AFf) is the free energy re- 
quired to form the defect system from a host (defect-free) 
lattice H . For a defect system X consisting of dopant 
atoms and associated lattice vacancies, which are formed 
by the removal of atoms of chemical species i, AFf 
becomes (for any temperature T within the harmonic 
approximation) : 



AFf = F (X) — 



cat. 



(1) 

Here F (X) is the total free energy for the defected sys- 
tem, and F (H) the total free energy of the defect-free 
host lattice. The chemical potentials of the dopant atom 
(D) and the atom(s) of chemical species i removed to 
create a vacancy in the defect lattice are [Ib and jii . The 
number of the dopant atoms inserted is no', n.i is the 
number of lattice atoms removed to create vacancies. Fi- 
nally, cat corresponds to charge correction terms which 
are added to equation [T] if the atoms removed or added 
have a charg o 29 ' 40 . 

The free-energies are defined as 

F(X) = E (X) + F vih (X) - TS conRg , (2) 

F{H) = E (H) + F vih (H) , (3) 

where E (X) is the total (DFT) energy for the relaxed 
defect system, and E (H) is the total (DFT) energy of 
the relaxed defect-free host lattice, F v ;b is the vibrational 
free energy contribution from the system X or host H as 
given by 
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where g(ui) is the (harmonic approximation) phonon den- 
sity of states for system X or H . It is important to note 
that here g(u>) has been normalized such that 
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Here Np is the number of phonon branches in the super- 
cell. As written, _F v ib also contains the zero-point energy 
Ezp = J Q g (oj) \hwduj contribution of the X or H sys- 
tem, viz., 
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The configurational entropy contribution S con fi g is given 
by S'config = kg In(fix)) where fix enumerates the pos- 
sible orientations of the defect (X) within the supercell. 
In both cases, our defects require the Xe atom to be on 



a lattice site (initially), and have a distinct orientation 
towards one of the nearest-neighbor lattice sites. This 
evaluates as: 



n 



X 



nin = 864, 



(7) 



where ni is the number of lattice sites in the supercell 
(216) and n a is the number of distinct orientations the 
Xe defect may have per lattice site (or number of nearest- 
neighbor C atoms, 4). While this works for Xe-3V, for 
Xe-V we must, as the final structure is significantly dif- 
ferent from the initial (see Sec. IIV|) . include a term to 
avoid double-counting: 



„ nin n 
ttxo-v = -77^ = 432. 



(8) 



The chemical potential of the carbon atoms removed 
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E DFT (H) + F vih (H) 
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where n is the number of C atoms in the defect-free su- 
percell of diamond. For this work we consider the forma- 
tion of a neutral dopant (Xe atom) in a covalent system 
(diamond), thus charge correction terms may be ignored. 

As we are considering two defected systems (Xe-V and 
Xe-3V) , in determining the relative thermodynamic sta- 
bility of the two defect configurations we need only con- 
sider the relative defect formation free energy 
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AFf (Xe-V) - AFf (Xe-3V) (10) 



allowing us to divorce common bulk properties from rel- 
ative stability calculations: 



AF 



REL 



F (Xe-V) - [F (Xe-3V) + 2/x c ] 



(11) 



As defined, AFf EL > means that the Xe-3V case is 
more stable than the Xe-V case. 

Note that the fixe term does not appear in the above 
expression, due to both our defects being comprised of 
one Xe dopant atom, and the /ixe terms cancelling. 



III. METHODOLOGY 

All DFT calculations were performed using the ah 
initio total-energy and molecular dynamics package 
VASP (Vienna ab initio simulation package) developed 
at the Universitat Wien 4 ^. We used the generalized- 
gradient approximation (GGA), with the PBE exchange- 
correlation functional 42 , and Plane Augmented Wave 
(PAW)— potentials for C and Xe. The plane-wave cut- 
off was set at 400 eV for all calculations. A 4x4x4 
Monkhorst-Pack fc-point mesh was used as this has been 
found to give accurate representation of the geometry and 
electronic structur o 29 ' 34 i 44 ' 45 . The bulk diamond lattice 
constant was checked using these parameters, found to 



Bond Length 


Relaxed 


Change from 


/ Angle 


amount 


initial 


Ci-Xe 


2.15A 


+39.6% 


ZCi-3-Xe-Ci- 3 


84.436° 


-22.9% 


ZC4-6~Xe— C4-6 


84.316° 


+40.5% 


ZCi-3-Xe-C4-6 


95.640° 


+ 6.27% 



TABLE I: Bond lengths and angles for relaxed 
split-vacancy Xe 



be 3.562 A, and this value was used for all subsequent 
calculations. 

All calculations were performed by constructing a 
3x3x3 simple cubic supercell of diamond (216 atoms), 
removing the appropriate number of carbon atoms to cre- 
ate vacancies and inserting Xe into a substitutional lat- 
tice site. The vacancy sites were always in first nearest 
neighbor positions to the Xe atom. For the defect dia- 
mond lattices, the Xe- vacancy defects were constructed 
such that they had trigonal crystallographic point group 
symmetry in accordance with experimental data31. 

In performing the geometry optimization the relax- 
ation was constrained so that the supercell remained cu- 
bic in shape. After relaxation the force on a C atom 
was typically less than 0.001 eV/A and in all cases the 
trigonal symmetry of the defect remained. 



IV. DEFECT GEOMETRIC STRUCTURE 
A. Split-vacancy interstitial Xe 

Upon relaxation the initial Xe- vacancy geometry (Xe s - 
V in Fig. [Ta|) . with all atoms occupying C lattice sites, 
formed a split-vacancy interstitial structure with the site 
symmetry of the Xe atom being C% v (Fig. Ilbl) . The Xe 
atom was then found to be at the centre of the divacancy, 
where it has D%d symmetry (still within the constraint of 
trigonal symmetry). In the relaxed structure, the Xe 
atom is six-fold coordinated with the nearest-neighbour 
C atoms, with bond lengths of 2.15 A. The bonds form 
two distinct sets of three (Ci_3 and C4_6); bond angles 
within sets are different from those between sets. More 
details may be seen in Table HI This structure is commen- 
surate with the Si— V defect in diamond studied by Goss 
et alM 



B. Three- vacancy defect (Xe s ~ 3V) 

The three- vacancy case (initial geometry shown in Fig. 
I2"a)) was also considered. Initially, the Xe is left with one 
nearest- neighbor C atom (Ci), and three directly adja- 
cent vacant sites. Of the surrounding atoms, nine (C2-10) 
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FIG. 1: (Color online) Split-vacancy Xe structure, 
viewed from [100] direction. C atoms shown in grey, 
nearest neighbors solid grey and numbered, Xe atom 
shown larger in blue, (a) Initial structure, (b) Relaxed 
structure. 



are also adjacent to the vacancies. Six of these (Cs_io) 
form a coplanar hexagon about the Xe atom. 

The relaxed geometry (Fig. I2b[) shows that the Xe 
atom migrates 0.38 A (out of the hexagon plane) to an 
interstitial position between the four vacant lattice sites 
(one being the initial position of the Xe atom) with the 
site symmetry of the Xe atom being C^ v . Ci_4 move 
outward by 0.22 A. The hexagon distorts slightly, with 
pairs of atoms expanding outwards (along their perpen- 
dicular bisector) by 8.66 A, maintaining symmetry. 
More details are located in Table |TTJ 
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FIG. 2: (Color online) Xe s -3V, viewed from 
approximately the [110] direction. C atoms shown in 
grey, nearest neighbors solid grey and numbered, 
6- member C ring shown in pink and numbered, Xe 
atom shown larger in blue, vacant sites shown as blue 
wire ellipse frames, (a) Initial structure, (b) Relaxed 
structure. 



V. DEFECT FORMATION ENERGIES 

Using the equations detailed in Section [Til the rela- 
tive defect formation free energy (AF^ EL from Eq. ITT]) 
was calculated over temperatures ranging from 0-1500K 
and can be seen in Fig. [3J The split-vacancy defect is 
2.59 eV more stable than the Xe-3V defect at all tem- 
peratures considered. The variation with temperature is 
small, about 3% of the total value, indicating that the 
difference in their behaviour with respect to temperature 
is small compared to the other terms in AFj lEL and that 
for the purpose of determining the qualitative stability 
ordering temperature effects are irrelevant. The relative 
stability of the split-vacancy defect is slightly enhanced 
at high temperatures. Neglecting Ezp (which comes at 
significant computational expense compared to a simple 
DFT calculation), leads to a relative zero-temperature 
defect formation energy of -3.06 eV, overestimating the 
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initial 


Ci-Xe 


2.13A 


+38.3% 


C2-4 - Xe 


2.39A 


- 5.28% 


Cs-io-Xe 


2.63A 


+ 4.47% 


C 5 -C 6 (7-8,9-10) 


2.67A 


+ 5.86% 


C 6 -C 7 (8-9, 10-5) 


2.52A 


+ 0.11% 


^-KJl Atj ^5-10 


80.334° 


-10 7% 


ZC 5 -Xe-C 6 (7-Xe-8, 9-Xe-10) 


60.878° 


+ 1.46% 


ZC 6 -Xe-C 7 (8-Xe-9, 10-Xe-5) 


57.247° 


- 4.59% 


ZC2-4-Xe-C5_io 


67.828° 


+13.0% 


^C2-4 _ Xe— C2-4 


63.129° 


+ 5.23% 



TABLE II: Bond lengths and angles for relaxed Xe-3V 
stability by wl8%. Finally, in Fig. @] we show the vi- 
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FIG. 3: Relative defect free energy of formation; 
F (Xc-V) - [F ( Xe-3V ) + 2/mo]. (Color online) 

brational free energy contribution to the total energy as 
a function of temperature for each defect. The Xe-3V 
case (dashed line) is vibrationally more stable at every 
temperature studied. At OK, the difference is 0.83 eV, 
and at 1500K, 0.25 eV. We may therefore infer that the 
energetic expense of removing two C atoms from the di- 
amond lattice outweighs the steric benefit of having the 
extra space for the Xe atom and unpaired C electrons to 
inhabit. 



VI. VIBRATIONAL DOS 

The vibrational properties of the systems were then 
calculated using lattice dynamics under the harmonic ap- 
proximation. The derivatives of the forces required to 
build the dynamical matrix were obtained using a finite- 
difference scheme in which individual nuclei were dis- 
placed by 0.04 A from their equilibrium positions and the 
resulting forces on all the nuclei calculated. Symmetry 
was used to reduce the number of displacements required 
to build the full dynamical matrix. The calculated forces 



5 



> 

0) 



40- 



20- 



0- 



200 400 



600 800 1000 
Temperature (K) 



1200 1400 



FIG. 4: Vibrational contribution to the free energy. 
Split- vacancy Xe defect (solid line), Xe-3V defect 
(dashed line). (Color online) 



were used in the program phon^ to construct and di- 
agonalise the dynamical matrix to obtain the vibrational 
Density of States (vDOS) and the thermodynamic prop- 
erties of the systems. 

Figure [5a] shows the total vDOS of bulk diamond 
(shaded) and for the Xe split-vacancy (black), averaged 
over the entire Brillouin zone. The region where the 
contribution from the Xe atom, the Xe-projected vDOS 
(dashed red) is significant is shown in Fig. 5b. Three 
strong Xe-localized modes are extant at the T point (dot- 
dashed green). As the site symmetry of the Xe atom 
and the crystallographic point group symmetry of the 
relaxed lattice are identical (C3,,), the vibrational modes 
involving the Xe atom will be of A- (non-degenerate) or 
.E-type (doubly-degenerate) symmetry. We find a non- 
degenerate mode with Ai-character at 218 cm" 1 whose 
eigenvector shows the displacement of the Xe atom to 
be in the [111] direction. This displacement compresses 
three of the C-Xe bonds (e.g. Ci_3) and stretches the 
other three. A doubly-degenerate mode of ^-character 
was found at 251 cm" 1 which is more prominent in the 
full vDOS than the A\ mode. All three modes are mu- 
tually orthogonal (see Fig. [5]), and are both IR- and 
Raman-active. For all three modes the relative displace- 
ment of any C atom is much smaller (by an order of mag- 
nitude) than the displacement of the Xe atom, showing 
that the vibrations are indeed almost entirely localized 
on the Xe atom. 

The corresponding total vDOS of the Xe-3V defect 
(black) and for bulk diamond (shaded) is shown in Fig. 
[7al (again averaged across the Brillouin zone). Again, 
the region where the Xe-projected DOS (dashed red) is 
significant is displayed in Fig. 7b. At the T point we find 
a non-degenerate mode at 245cm _1 (dot-dashed green) 
which displaces the Xe atom along the [111] direction 
(along the Xe-Ci bond). The relative displacement of 
any C atom is again smaller (by an order of magnitude) 
than the displacement of the Xe atom. As this mode is 
non-degenerate, it corresponds to an Ai-type mode which 
is both IR- and Raman-active. 

From 400 cm -1 , both vDOS spectra are largely un- 
changed from the bulk, apart from some softening of 
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FIG. 5: (Color online) Vibrational density of states for 
the relaxed split-vacancy Xe defect. Total vDOS (solid 

black line), Xe-projected vDOS (dashed red line), 
Gaussian-smeared Xe-projected vDOS at T (dot-dashed 
green line), bulk diamond vDOS (shaded background), 
(a) Full calculated vDOS. (b) Region where 
Xe-projected vDOS is significant. 



peaks and a slight shift from about 900 cm -1 onwards. In 
particular, the sharp peak near 1300 cm -1 is significantly 
lower. 

The predicted Raman signals should be visible and dif- 
ferentiable at sufficient doping density, especially if ap- 
propriate line-fitting models are used^i. Some indications 
of them are present in the PL phonon sidebands seen in 
Martinovich et al^, but a clearer signal could be ob- 
tained by looking at the sidebands nearest the pump laser 
frequency. This may reduce any frequency-shift effects of 
the electronic level excitation/de-excitation. Also, if the 
center in question happens to be the optically active one, 
pumping its optical emission could set up a dipole which 
may enhance the Raman signal, making detection easier. 
Unfortunately, none of the currently published experi- 
mental work on Xe centers provide spectra sufficiently 
close to the pump laser frequency to compare with-^"— . 



VII. CONCLUSIONS 

We have presented the results of Density Functional 
Theory simulations of Xe in diamond and draw several 
conclusions. First, the split-vacancy neutral defect is the 
most stable of those considered, at all temperatures. Sec- 
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FIG. 6: Displacements of the Xe atom under the two 
strongly Xe-localized vibrational modes in the 
split- vacancy Xe defect. Aj mode (red arrow), 
doubly-degenerate E mode (green arrow). (Color online) 



ond, the cost of removing two extra C atoms from the lat- 
tice outweighs the vibrational benefit of the extra space. 
Finally, both defects have IR- or Raman-active vibra- 
tional modes at T which may aid in further identification 
of the structure of the defects studied by Bergman et 
alM&. 
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FIG. 7: (Color online) Vibrational density of states for 
the relaxed Xe-3V defect. Total vDOS (solid black 
line), Xe- projected vDOS (dashed red line), 
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